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We analyze the competition between antiferromagnetism and superconductivity 
in the two-dimensional Hubbard model by combining a functional renormalization 
group flow with a mean-field theory for spontaneous symmetry breaking. Effective 
interactions are computed by integrating out states above a scale Amf in one-loop 
approximation, which captures in particular the generation of an attraction in the d- 
wave Cooper channel from fluctuations in the particle-hole channel. These effective 
interactions are then used as an input for a mean-field treatment of the remaining 
low-energy states, with antiferromagnetism, singlet superconductivity and triplet vr- 
pairing as the possible order parameters. Antiferromagnetism and superconductivity 
suppress each other, leaving only a small region in parameter space where both 
orders can coexist with a sizable order parameter for each. Triplet vr-pairing appears 
generically in the coexistence region, but its feedback on the other order parameters 
is very small. 



Soon after the discovery of high-temperature superconductivity in cuprate compounds, 
Andersor>i pointed out that the essential physics of the electrons in the copper-oxide 
planes of these materials could be described by the two-dimensional Hubbard model. The 
model describes tight-binding electrons with a local repulsion ?7 > 0, as specified by the 
Hamiltonian 



in standard second quantization notation. A hopping amplitude —t between nearest 
neighbors and an amplitude —t' between next-nearest neighbors on a square lattice leads 
to the dispersion relation ek = — 2t (cos kx + cos ky) — At' cos kx cos ky for single-particle 
states. 



PACS: 71.10.Fd, 74.20.-z, 75.10.-b 



I. 



INTRODUCTION 
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In agreement with the generic phase diagram of the cuprates, the Hubbard model is 
an antiferromagnetic insulator at half-filling (provided t' is not too big), and is expected 
to become a d-wave superconductor away from half-filling in two dimensions already for 
quite some time.- In particular, the exchange of antiferromagnetic spin fiuctuations has 
been proposed as a plausible mechanism leading to d-wave pairing.-i^i^ It turned out 
to be very hard to detect superconductivity in the Hubbard model by exact numerical 
computation,^'^ as a consequence of finite size and/or temperature limitations. 

The tendency toward antiferromagnetism and d-wave pairing in the two-dimensional 
Hubbard model is present already at weak coupling. However, conventional perturbation 
theory breaks down in the most interesting density regime, since competing infrared diver- 
gences appear as a consequence of Fermi surface nesting and van Hove singularities.-'^'^ 
A controlled treatment of these divergences can be achieved by a renormalization group 
(RG) analysis, which takes into account the particle-particle and particle-hole channels 
on equal footing. 

A suitable framework for a systematic RG analysis of the 2D Hubbard model is pro- 
vided by the so-called exact or functional RG.— In this approach, fermionic fields in a 
functional integral representation of the model are integrated successively by descend- 
ing step by step in energy scales. This can be formulated as an exact hierarchy of fiow 
equations for the effective interactions. The energy scale of the fields. A, is the fiow 
parameter. Several groups have computed the flow of effective two-particle interactions 
for the two-dimensional Hubbard model, using various versions of the functional RG 
in one-loop approximation.— '^'^'^ Antiferromagnetic and superconducting instabilities 
were detected from the flow of the corresponding susceptibilities. 

At sufficiently low temperatures, and in particular at T = 0, the effective two-particle 
interaction obtained from a one-loop approximation diverges at a flnite scale Ac > 0, 
that is, before all flelds have been integrated out. Hence, one is running into a strong 
coupling problem in the low-energy limit, even in the case of a weak bare interaction. 

If the vertex function diverges only in the Cooper channel, driven by the particle- 
particle contribution to the flow, the strong coupling problem emerging in the low-energy 
region can be controlled by exploiting Ac as a small parameter.— The scale Ac is exponen- 
tially small for a small bare interaction. The formation of a superconducting ground state 
can then be described essentially by a BCS mean-field theory with renormalized input 
parameters. In the two-dimensional Hubbard model, an instability where the interactions 
diverge only in the Cooper channel is realized at sufficiently small [/, if the Fermi surface 
stays away from van Hove points.— 'i^'^'i^ At van Hove filling the effective interactions di- 
verge also in other channels even in the weak coupling limit. Although superconductivity 
then has to compete with other instabilities, the one-loop flow indicates that it is still the 
leading instability for a moderate t' 7^ and sufficiently small U .—'^^^ 

In principle, spontaneous symmetry breaking can be handled within the functional RG 
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framework by adding an infinitesimal symmetry breaking term at tlie beginning of the 
flow, which is then promoted to a finite order parameter at the scale Ac.— So far, this 
approach has been worked out in practice only for mean-field models. Order parameter 
fluctuations are most conveniently treated by introducing appropriate bosonic fields, as 
discussed recently for antiferromagnetic order in the half-filled Hubbard model.— 

In case of competing order parameters, such as antiferromagnetism and d-wave super- 
conductivity, a full RG treatment of spontaneous symmetry breaking and order parameter 
fluctuations is a rather ambitious long-term goal. In the present work, we explore a simpler 
alternative and combine the RG with a mean-field (MF) theory of symmetry breaking. 
In this RG+MF approach the one-loop flow is stopped at a scale Amf that is small com- 
pared to the band width, but still safely above the scale Ac where the two-particle vertex 
diverges. At this point the vertex has developed already a pronounced momentum depen- 
dence, reflecting in particular antiferromagnetic and superconducting correlations. The 
integration over the remaining modes, below Amf, is treated in a mean- field approxima- 
tion allowing, in particular, antiferromagnetic and superconducting order. Low-energy 
fluctuations are thereby neglected. The mean-field Hamiltonian is defined on a restricted 
momentum region near the Fermi surface, with |ek — /i| < Amf, and the effective interac- 
tions entering the mean-field equations are extracted from F'^^^'. A very short account of 
some RG+MF results has appeared recently in Ref. 0. 

Our theory extends previous mean-field treatments of antiferromagnetism and d-wave 
superconductivity in two-dimensional Hubbard-type and t-J models, where the effec- 
tive interactions were specified by an ad hoc ansatz or identified with bare microscopic 
interactions,— i^^i^^i^^ while we compute the effective interactions by integrating out fluc- 
tuations. 

In Sec. II we will review the functional RG and the structure of the flow equations for 
the effective interactions on one-loop level, focussing on the Wick ordered version used in 
the present work. The mean-field equations for antiferromagnetism and superconductivity, 
including a possible coexistence of both, will be derived in Sec. III. Results from the 
combined RG+MF analysis for the 2D Hubbard model will be presented in Sec. IV. In 
particular, we will analyze which phases are stabilized for various choices of hopping, 
interaction and density, and we will show results for the size and shape (momentum 
dependence) of the relevant order parameters. We conclude with a summary of the main 
results in Sec. V. 

II. RENORMALIZATION GROUP 

All versions of the functional RG for interacting Fermi systems are variants of 
Wilson's^^ momentum shell RG, where fermionic fields are integrated out by descending 
successively from the modes with highest energy down to the Fermi surface.^'' This type 
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of RG is also the basis for important rigorous work on two-dimensional Fermi systems.— 
The successive integration of modes can be formulated as an exact hierarchy of flow equa- 
tions for effective interactions (one-particle, two-particle etc.).— A similar hierarchy 
can be obtained by using the interaction strength instead of a cutoff as flow parameter.— 
We focus on the Wick-ordered version of the functional RG, which we will use in our 
calculations. 



A. Wick ordered flow equations 



Before turning to approximations, let us first sketch the structure and derivation of 
the exact flow equations (for details, see Salmhofer— and Ref. Il3l ). The flow parameter A 
is introduced as an infrared cutoff for the bare propagator, such that contributions from 
states with momenta obeying |ek — /i| < A are suppressed. All Green functions of the 
interacting system will then depend on A, and the true theory is recovered only in the limit 
A ^ 0. The RG equations are most conveniently obtained from the effective potential V^, 
which is the generating functional for connected Green functions with bare propagators 
amputated from the external legs. Taking a A-derivative one obtains an exact functional 
flow equation for this quantity. Expanding on both sides of the flow equation in powers 
of the fermionic fields (i.e. Grassmann variables), and comparing coefficients, one obtains 
the so-called Polchinski equations^ for the effective m-body interactions. These equations 
were used by Zanchi and Schula^^ in their RG analysis of the two-dimensional Hubbard 
model. An alternative expansion in terms of Wick ordered monomials of fermion fields 
yields flow equations for the corresponding m-body interactions V.^ with a particularly 
convenient structure (see Fig. 1).^- The flow of is given as a bilinear form of other 



n, i Perni. 





FIG. 1: Diagrammatic representation of the flow equation for Vj^ in the Wick ordered version 
of the functional RG. The line with a slash corresponds to dG^^ /dA, the others to Gq^; all 
possible pairings leaving m incoming and m outgoing external legs have to be summed. 



n-body interactions (at the same scale A), which are connected by lines corresponding to 
the propagator 

.<A,. u^ ©(A-I^kl) 



(2) 



iko - 

where = — and one line corresponding to SaGq "^(^o, k). Note that the support 
of G'g^(A;o;k) is restricted to momenta with |^k| below the scale A. These soft mode 
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propagators come into play via the Wick ordering, which is given by (scale dependent) 
contractions with Go^(A;o,k)- For A 0, the effective interactions remain unaffected 
by the Wick ordering, since GQ^{k) vanishes in that limit. For small A, the momentum 
integrals on the right hand side of the flow equation are restricted to momenta close to 
the Fermi surface (see Fig. 2). 
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FIG. 2: (Color online) A typical Fermi surface (bold line) and the support of the propagator 
GQ^{ko,'k) (dotted region) in the first quadrant of the Brillouin zone. The dashed lines mark 
the boundaries of the "patches" used for the discretization of the momentum dependence of the 
effective two-particle interaction. 

With the initial condition V"^" = bare interaction, where Aq = max|^k|; the above 
flow equations determine the exact flow of the effective interactions as A sweeps over the 
entire energy range from the band edges down to the Fermi level. Since the Wick ordered 
flow below scale A involves only low-energy states with |^k| < A, it yields a continuous 
sequence of effective low-energy models with effective interactions acting on a restricted 
momentum space. 



B. One-loop flow 

To detect instabilities of the system in the weak-coupling limit, it is sufficient to trun- 
cate the infinite hierarchy of flow equations described by Fig. 1 at second order in the 
effective two-particle interaction V.^. Contributions to the flow of V.^ involving effective 
three-particle interactions and higher m-body terms are at least of third order in V^^. The 
leading contribution from the effective one-particle interaction can be absorbed by a 
simple shift of the chemical potential. The remaining influence of on the flow of 
is of third order an can thus be neglected. 

The flow of the effective two-particle interaction V-^ is thus given by one-loop terms 
involving only V^, and no other m-body terms. In the following, we write F^ instead of 
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V.^ for the two-particle interaction. Putting arrows on the hnes to distinguish creation and 
annihilation operators, one obtains the diagrammatic representation of the flow equation 
shown in Fig. 3. To write the momentum, energy and spin dependences of F"^, we collect 
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FIG. 3: Flow equation for the effective two-particle interaction F^ in one- loop approximation 
with the particle-particle channel (PP) and the two particle- hole channels (PH and PH'). 

k = (/cq, k) and a in a single variable K. In T^{K[, K2, Ki, K2), the variables Ki and K2 
refer to incoming, the variables K[ and K2 to outgoing particles. With this notation, the 
explicit one-loop flow equation for F'^ reads^ 

d 



dA 



d 



^i:,aA[Go-(A-)G„-(A-' 



K,K' 



X 



(3) 



[,f\K[, K'2, K, K') T\K, K'- K,, K^] 
- V\K[, K; K,, K') V\K' , K'^, K, K^) 
+ V\K'2, K- K,, K') V\K\ K[- K, K^) 

where L is the number of lattice sites, and (3 the inverse temperature. The three terms on 
the right hand side are the contributions from the particle-particle channel (PP) and the 
two particle-hole channels (PH and PH'). For translation invariant systems momentum 
conservation implies that V^[K[, K'2] Ki, K2) 7^ only if ki + k2 = k[ + k'2, so that the 
sum over k and k' in ([3]) is reduced to a single energy-momentum sum. 

For a spin-rotation invariant system the spin structure of the two-particle interaction 
can be written as 



r\K[, K'2; K,, K2) = F^(A;;, k'2; ku ^2) ^.j,.-.,,., + Pf (fc'i, k'2; k,, ^2) T^^.-o 
where 



rp 



(4) 



(5) 
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are the projection operators on singlet and triplet states in a two-particle spin space, 
respectively. The antisymmetry of with respect to Ki <-> K2 or K[ ^ K'2 implies that 

is symmetric and antisymmetric under exchange of the variables ki and k2 or k[ 
and /cg- Inserting the above decomposition in singlet and triplet terms, the spin sums can 
be easily carried out and one obtains a coupled set of equations for and F^. 

It is clearly impossible to solve the flow equations with the full energy and momentum 
dependence of the two-particle interaction, since F^ has three independent energy and 
momentum variables. The problem can however be simplified by ignoring dependences 
which are irrelevant (in the RG sense) in the low energy limit, namely the energy depen- 
dence and the momentum dependence normal to the Fermi surface.— Hence, we compute 
the flow of the two-particle interaction at zero energy and with at least three momenta on 
the Fermi surface (the fourth being determined by momentum conservation). On the right 
hand side of the flow equation we approximate the interaction by its zero energy value 
with three momenta projected on the Fermi surface (if not already there). Note that this 
projection procedure is exact for the initial two-particle interaction, since the Hubbard 
interaction is momentum and energy independent. The remaining tangential momentum 
dependence is discretized. The momentum-dependence of the two-particle vertex is thus 
approximated by a step function which is constant on "patches" (sectors) in the Brillouin 
zone.— '^'^ Here we use patches deflned by straight lines connecting the points (0, 0) and 
(tt, tt) to the magnetic Brillouin zone boundary (also known as "umklapp surface"), see 
Fig. 2. 

Neglecting the energy dependence of F^, the Matsubara sum in the flow equation can 
be done analytically. Due to the sharp momentum cutoff the momentum integral can be 
easily reduced to a one-dimensional integral. The latter, and the integration of the flow, 
has to be performed numerically. 

We will stop the flow at a scale Amf well above the scale Ac at which F^ diverges. 
The remaining degrees of freedom will be treated in a mean-fleld approximation to be 
described in the next section. The output of the flow, F^^i", is the input interaction for 
the mean-fleld theory. 

III. MEAN-FIELD THEORY 

The divergence of the effective two-particle interaction F^ as a function of decreasing A 
signals an instability leading to an ordered phase via spontaneous symmetry breaking.— 
The simplest way to treat spontaneous symmetry breaking is provided by a mean-fleld 
approximation. For a speciflc type of order, only interactions with a very restricted choice 
of momenta are picked up by the mean-fleld theory. For example, in the BCS mean-fleld 
theory of superconductivity, the only relevant interaction processes involve particles with 
vanishing total momentum (ki -|- k2 = 0). Vice versa, in the thermodynamic limit the 
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mean-field approximation provides the exact solution of a reduced Hamiltonian, in which 
the interaction is restricted to the relevant momenta from the beginning. 

A. Reduced Hamiltonian 

Here we specify our reduced Hamiltonian, and link the interaction terms to the effec- 
tive two-particle interaction obtained from the RG. In addition to interactions driving 
antiferromagnetism and superconductivity, we will also include several other interaction 
terms, since they do not lead to major complications, make the general structure of the 
theory more transparent, and may be useful for a more general analysis in the future. 

The reduced Hamiltonian has the form 

= Ho + Hl"^ , (6) 

where Hq = ^ n^^ is the kinetic energy and H^'^'^ the reduced interaction. The latter 
contains four terms, 

which we now specify one by one. 

The first term is a density-density interaction with zero momentum transfer (forward 
scattering) , 

k,k' cr,cr' 

This term captures Fermi liquid interaction effects, and can lead to charge or spin den- 
sity instabilities such as phase separation, ferromagnetism, or a d-wave Pomeranchuk 
instabilityi^*^. 

The second term is a density-density interaction with momentum transfer Q = (tt, tt), 

^r'^ = ^EEf^kk'<<^<v', (9) 

k,k' cr,cr' 

with n^o- = ojccr'^k+Q.cr- Note that k + Q = k — Q. This term drives charge or spin density 
wave instabilities with a wave vector Q and arbitrary form factors (s-wave, d-wave etc.), 
including, in particular, antiferromagnetic order. 

The third term is a singlet pairing interaction between particles with total momentum 
zero (Cooper channel), 

k,k' 

with = a_kj,ctkT ^-^id pj^ its hermitian conjugate. There is no factor 1/2 here, since spin 
variables are fixed. This term drives spin-singlet superconductivity. 
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Finally, we include a triplet pairing interaction between particles with a total momen- 
tum Q, 

Hr = jJ2^^^'pVpi'^ (11) 

k,k' 

with p'^ = aQ_k,j,akT and p^^ its hermitian conjugate. In case of coexistence of antifer- 
romagnetism and superconducitivity, this term leads to a condensate of pairs with total 
momentum QM^^ 

We use the mean-field theory as an approximate treatment of the fermionic degrees 
of freedom below the energy scale Amf- Hence, all momenta in the above terms are 
restricted to a shell around the Fermi surface given by |ek — /u| < Amf- For the dispersion 
relation in Hq, we simply use the bare one, that is, we do not keep track of self-energy 
contributions which may renormalize ek- Self-energy corrections beyond those which can 
be absorbed in a shift of /i appear only at second order in the interactions, and may 
therefore be neglected at weak coupling. 

The coupling functions in Hl'^'^ are extracted from the effective two-particle interaction 
pAMF Denoting the static (frequency independent) two-particle interaction with generic 
momenta by T^^^ (k[a[, ^(72; kicii, k2cr2), we have 

^kk' = I r^^^ka, kV; ka, kV) , (12) 
f^kk' = I r^^^ka, kV; k + Q a, k' - Q ^0 , (13) 
lick' = ^r^-(k,-k;k',-k'), (14) 

^kk' = I r^-(k, Q - k; k', Q - k') . (15) 

All these coupling functions are real-valued. Note that there is an overlap between 
the different coupling functions above for special choices of momenta. For example, for 
ki = k'^ = — k2 = — k2 the vertex belongs to the Cooper and forward scattering chan- 
nels simultaneously. However, these special momentum sets have zero measure in the 
thermodynamic limit. 
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B. Mean-field decoupling 

The mean-field theory involves the following mean fields 

/^k. = ^J^i^kk'Kv) (16) 

k'cr' 
k'cr' 

Ak = -^5^Fkk'(Pk') (18) 
k' 

= ^E^kk'(Pk')- (19) 

k' 

The density field Dko- is real. The density field D^^ can be complex and obeys the relation 
DZ+ci,a = Di:, which follows directly from {K_^^^J = (n^J* and t/^^Q,k'+Q = U^i^, . The 
pairing fields Ak and A'^ are generally complex. 

The operator products in the interaction terms have the form b'b, where b and b' are 
products of two elementary Fermi operators. In the mean-field decoupling, one approx- 
imates b'b ^ {b') b + (6) b' — {b') (b). The fluctuation term {b' — {b')) (b — (b)) is thereby 
neglected. This yields 



Dka ^ka - ^ E Dka (^ka) , (20) 
k.cr k,cr 

E^k.<.-^E^^^>L), (21) 

k,o- k,cr 

Hf'' ^ E ^kPk + 5^ AkpL - J] A* (pk) , (22) 

k k k 

Hr ^ E^k>k+E^kPk^-E^r(Pk)- (23) 

k k k 

For the above reduced interactions this approximation is exact in the thermodynamic 
limit. This can be seen quite easily by estimating the size of Feynman diagrams with a 
reduced interaction or, alternatively, by using a Hubbard Stratonovich decoupling in a 
functional integral representation.— 

The Hamiltonian is now quadratic in the fermions and can therefore be diagonalized. 
We define = — fiN, where is the reduced Hamiltonian in mean-field 

approximation, and = ^ nko- is the total particle number operator. Introducing the 
Nambu spinor 

ak = (flkT, a^ki' ^k+Q,T5 «lk-Q,i) ' (24) 



11 



Ml 



(26) 



the mean-field Hamiltonian can be written as 

= J2slM^Si^ + K\ (25) 

k 

where K"^ collects all c-number terms, that is, terms without operators, and the quadratic 
form with the matrix TVlk all quadratic terms. Here the momentum sum extends only 
over momenta in the magnetic Brillouin zone. The matrix A^k is given by 

^k -^-k - ^-ki ^r+Q -D-k-Q,i 

-Dk+Q,T ^k+Q ^k+Q + ^k+Q,T ^k+Q 

V Ak* -^-ki Ak+Q -e-k-Q - /^-k-Qa / 

Note that Aik is hermitean since -Dk+Qcr — ^ka- The momenta of the mean fields in 
Aik are restricted by the cutoff: Dko- and Ak are restricted by the condition |^k| < ^mf, 
while and AJJ contribute only if the two conditions |^k| < ^mf and l^k+ql < ^mf are 
satisfied simultaneously. The matrix TWk has a 2 x 2 block structure 

where and Ai^ are 2x2 matrices. 

The c-number term, K'^ contains a contribution from commutators, ^k(^k + -Dki), 
in addition to the c-number terms originating from the mean-field decoupling of the 
interaction terms. 

The matrix A^k can be diagonalized by a unitary transformation U^: 



alMk ak = XI 4-^k ak , (2^ 



k k 

where ak = Wkak, a|. = a|^W,[, and Ai^ = UkJ^k.Ul is a diagonal matrix 

Mk = dmg{El,ElElEt) (29) 

with the eigenvalues E^, A = 1, 2, 3, 4, as entries. The transformed creation and annihila- 
tion operators obey fermionic anticommutation relations, {a^,a^,^} = 5xx' S^k/ , since 
is unitary. The block structure ( 1271) implies that the matrices A^k and Alk+Q have the 
same eigenvalues, which can thus be labelled such that -^k+q — -^k- 

C. General gap equations 

The selfconsistency or gap equations can be derived by minimizing the grand canonical 
potential 

n = K" -Tj2^n(^l + e-^^^^ (30) 

k,A 
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with respect to the expectation values (nko-), {n^^^), (pk), {p'Q, or complex conjugates. 
As shown in Appendix A, the gap equations can be expressed in terms of minors of the 
matrices 

Mi, = M^-E^l. (31) 

The minor det(Mj^'-'"' ) is the determinant of the matrix obtained from by deleting 
the j-th row and the j'-th column. Minimizing with respect to (nko-) yields 



1 V- V- ^kk' det (M^;") - F^'i^, det {M^!^') ^ . ^ 1 .1 



^i^V E,det(M,V') 



(32) 



k' 



where / is the Fermi function. The last term in (132!) stems from the contribution 
to K'^ mentioned above. Minimizing with respect to (T^ko-) yields 

E,det(M,y) 

Minimization with respect to (pk)* yields 



A^ = 4E^-E^^^%xlT/(i^^0, (34) 
VE«det(Mj,, ) 

while minimizing with respect to (pk) would just yield the complex conjugate of the above 
equation. Minimization with respect to (p^)* yields 

AJ^-4E^-E^^^ST/(i^^0, (35) 
X Eidet(Mj^; ) 

and minimization with respect to {p'^) again the complex conjugate. Note that the mo- 
mentum sums in the above gap equations extend over the full (not just magnetic) Brillouin 
zone within the limits imposed by the cutoff. 

The denominator in the gap equations can be expressed in terms of the eigenvalues as 
(see Appendix A) 

5^det(M,^'") = n(^k -^k)- (36) 

I X'^X 

The electron density n = —L~^dQ/dfi can be written as 

, l^^ det(M,V^)-det(M,y^) 
^T^x E«det(Mk; ) 

The derivation of this relation is almost identical to the one of the gap equation for -Dkcr- 
For a numerical evaluation it is sometimes convenient to rewrite the ratios of determi- 
nants in the above expressions in terms of the eigenvectors of A^k-— 
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D. Gap equations for antiferromagnetism and superconductivity 

We are mainly interested in the interplay of antiferromagnetism and superconductivity 
in the two-dimensional Hubbard model. In the numerical solution of the mean-field equa- 
tions we will therefore ignore -Dkcr. Furthermore, we will restrict D^^ to the spin density 
wave channel, that is, D^.^ = — -D^i = ^k- The mean-field Ay^ is related to the staggered 
magnetisation rrik = n'^-^ — n^^ by 



^k = ^5^f/ifk'(^k'), (3^ 
k' 



where U^y., = ^{U^^, — U^l^,'^). For a conventional commensurate antiferromagnet, is 
real. The matrix Ai]i then assumes the simpler structure 



Mv 



( a Ak Ak K \ 

A* -Ck A-Q Ak 

^k AJJ^Q ^k+Q Ak+Q 

V Ar Ak A*+Q -a+Q / 



(39) 



We have used reflection symmetry to replace ^_k by ^k etc. 

In principle, the coexistence of antiferromagnetic order and superconductivity leads 
automatically to a finite vr-pairing field.— i^iiM To see this, consider the right hand side of 
the gap equation ([35D for Aj^. The relevant minor in the numerator can be written as 

det(Mk"''^) = AkAk(^^ - a+q) + Ak+QAk(^^ + a) (40) 

in the limit Aj^ 0. This is nonzero if both Ak 7^ and Ak 7^ 0. In other words, the co- 
existence of antiferromagnetic order and superconductivity generates a finite expectation 
value (Pk)- For a nonzero interaction in the vr-pairing channel, V^k'' this leads to a finite 
TT-pairing field Aj^. Since both Ak and Ak are odd under spin flips, the spin flip symmetry 
of the expectation value (pj^) generated by Ak and AJ^ is even, corresponding to triplet 
pairing. This is why the interaction V^, was extracted from the triplet component of 
the vertex in Eq. (15). In case of coexistence of superconductivity with a charge (instead 
of spin) density wave, the generated vr-pairing would be odd under spin flips, and thus 
couples to the singlet vertex. Combined with the antisymmetry under particle exchange 
the spin flip symmetry of the vr-pairing field yields the relation AJ^ = — Aq_j^. 

At and near half-filling, the dominant ordering tendencies of the two-dimensional Hub- 
bard model are antiferromagnetism and singlet superconductivity. Under point group 
transformations the antiferromagnetic order parameter Ak has s-wave symmetry, while 
the superconducting gap Ak has d2:2_y2-symmetry. From Eq. (HOl) one can thus deduce 
that AJ^ also has d^2_j^2-symmetry. Note that the latter is not in conflict with the triplet 
spin symmetry of the vr-pair, since the total momentum of the pair is Q = (vr, vr), not 
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zero. The effective interaction T of the Hubbard model turns out to be very small in 
the triplet vr-pairing channel, and its d-wave component is repulsive. Hence the vr-pairing 
field AJ^ is also very small compared to the other order parameters, such that its feedback 
into the gap equations for and can be neglected. 

IV. RESULTS FOR THE TWO-DIMENSIONAL HUBBARD MODEL 

In this section we present results obtained from the renormalized mean-field theory on 
antiferromagnetic order and superconductivity in the repulsive two-dimensional Hubbard 
model. We will first discuss antiferromagnetic states in the absence of superconductivity, 
then d-wave superconductivity in the absence of antiferromagnetic order, and finally the 
full theory allowing for coexistence. Superconductivity in the attractive Hubbard model 
is discussed briefly in Appendix B. All results will be presented for t = 1, that is, all 
quantities are given in units of t. 

A. Antiferromagnetism 

In the absence of any other symmetry breaking, the gap equation for the antiferromag- 
netic order parameter can be written as 

= Z E ^kk' IfiK' - fiE^')] (41) 

k' k k 

with the two quasi-particle energy branches 

^± = §}i±^ ^k+Q)^ + Ai . (42) 

For a dispersion relation due to nearest and next-to-nearest neighbor hopping one has 

E^ = ^i±\/{^iy+Ai-^i, (43) 

where ej^. = — 2t(cos + cos ky) and = —At' cos k^ cos ky . 

The antiferromagnetic order parameter affects the quasi-particle energies most strongly 
near the magnetic Brillouin zone boundary, since vanishes there. The dispersion of the 
two branches along this line is shown in Fig. 4 for a constant (momentum indepen- 
dent) and t' < 0. Note that the momentum dependence of A]^ is indeed rather weak in 
our RG+MF calculations. Within plain mean-field theory applied to the Hubbard model 
one has U^^, = —U/2 so that is completely momentum independent. The antiferro- 
magnetic mean-field theory has been applied to the two-dimensional Hubbard model in 
several earlier works.— i2^»22i>^ 
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FIG. 4: (Color online) Dispersion of quasi-particle energies E^^ on the magnetic Brillouin zone 
boundary ("umklapp surface") in the antiferromagnetic state for a momentum independent Ai^ 
and t' < 0. 

The two branches and are separated by a global energy gap if |Ak| > 2|t'|. If the 
chemical potential lies in that gap, such that E^ > and E^ < for all k, the electron 
density is at half-filling and the system does not exhibit any gapless excitations. Gapless 
excitations do exist if E^ or E^ has zeros. The equations E^ = define effective Fermi 
surfaces of the antiferromagnetic state. From Figs. 4 and 5 one can see that for t' < 
and close to half-filling the effective Fermi surfaces enclose hole pockets in the branch E^ 
around (7r/2,7r/2) and electron pockets in E^ around (vr, 0) and (0,7r). Similar effective 




FIG. 5: (Color online) Effective Fermi surfaces in the antiferromagnetic state for t' = —0.2, 
/i = —0.6, and various momentum independent choices of the antiferromagnetic gap function 
Ai^: (solid line), 0.15 (dashed line), 0.3 (dashed-dotted line). 

Fermi surfaces are obtained in a d-density wave state.— For a momentum independent 
antiferromagnetic gap function the separatrices between states with hole pockets, 
electron pockets and fully gapped states are given by simple linear combinations of the 
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parameters t', /i and A. For a fixed t' < the various regimes in the plane spanned by fi 
and A are shown in Fig. 6. Note the special point given hj fii/2 = A = 2t' at which the 



hole-pocket 



0' ^ 



fully gapped 



electron-pocket' 



electron- and 
hole-pocket 



1^ 



FIG. 6: (Color online) Topology of the effective Fermi surfaces in the plane spanned by the 
chemical potential /i and the antiferromagnetic gap A. No Fermi surface exists in the fully 
gapped regime. 

separatrices for electron and hole pockets cross each other. For /i < /Xvh and a small A^ 
the effective Fermi surface is closed around (0, 0) and does not intersect with the magnetic 
Brillouin zone boundary. 

A typical mean-field result for the //-dependence of the antiferromagnetic gap in the 
Hubbard model is shown in Fig. 7. The coupling strength U has been chosen sufficiently 

0.5 



0.4 



a. 0.3 

us 

^ 0.2 



0.1 




-0.6 



-0.4 



-0.2 



FIG. 7: (Color online) Mean-field solution for the antiferromagnetic gap as a function of fj, for 
the Hubbard model with t' = —0.2 and U = 2.25. Regions with different Fermi surface topology 
are indicated as in Fig. [6l 



large to stabilize fully gapped half-filled solutions for /i near yUi/2- The antiferromagnetic 
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order parameter vanishes discontinuously at the edges of the antiferromagnetic region. 
The effective Fermi surface of the antiferromagnetic states consists (exclusively) of hole 
pockets for hole doping and of electron pockets for electron doping. For other finite values 
of t' the /i-dependence of is qualitatively the same for sufficiently large U. 

For t' = half-filled fully gapped antiferromagnetic solutions are stabilized for any 
(arbitrarily small) U, but no stable antiferromagnetic solutions away from half-filling 
exist. Solutions of the self-consistency equations leading to densities away from half- 
filling correspond to maxima (instead of minima) in the free energy and are therefore 
physically irrelevant.— If a certain density near half-filling is enforced, the system will 
thus phase separate in a half-filled antiferromagnetic and a non-half-filled paramagnetic 
region. 

We now turn to results obtained from the combined RG+MF theory described in the 
preceeding sections. The effective interaction f/^j^, driving the antiferromagnetic order is 
computed from the RG flow integrated down to a cutoff scale Amf, and the momentum 
space entering the mean-fleld equations is correspondingly restricted by this cutoff. Since 
the mode elimation is done only approximately both above and below Amf, the choice 
of Amf will affect the accuracy of the flnal results. In Fig. 8 we show the dependence 
of the antiferromagnetic gap function on Amf for the perfect nesting case (t' = at 
half-filling) with a relatively weak bare coupling strength U = 2. The different curves 
correspond to momenta k in six different patches interpolating between the or ky axis 
(patch 1) and the Brillouin zone diagonal (patch 6), see also Fig. 2. The patch-dependence 
is relatively weak, reflecting the weak momentum dependence of Ak. Fluctuations cap- 

1 



0.8 

0.6 



0.4 
0.2 



"01234 

FIG. 8: (Color online) Antiferromagnetic gap as a function of the cutoff Amf for the Hubbard 
model with pure nearest neighbor hopping {t' = 0) and [/ = 2 at half-filling. Due to the weak 
momentum dependence of the curves corresponding to different patches are very close to each 
other. The intersection of the largest ^k (patch 1) with the straight line given by Amf = 2^k 
yields at suitable choice of Amf leading to a reasonable estimate of the true gap size. 
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tured by the one-loop RG flow reduce the size of compared to the plain mean-field 
result (corresponding to Amf = Aq = 4). This reduction is expected and can also be 
obtained by a suitable perturbation theory for the symmetry-broken state.— Within 
mean-field theory, the effective interaction is driven exclusively by one-loop diagrams 
in the antiferromagnetic particle- hole channel, leading to an enhancement compared to 
the bare interaction. In the full one-loop flow other contributions, especially from the 
particle-particle channel, reduce this enhancement, leading to a smaller Ak. 

The size of Ai^ saturates at low Amf until it is pushed upwards by a strong divergence 
at a cutoff scale of the order of the gap size just before the upturn. This divergence 
reflects the breakdown of the one-loop approximation at the energy scale of symmetry 
breaking. The renormalized interaction becomes strong at that scale and the influence of 
the antiferromagnetic order parameter on the flow becomes crucial. The latter is captured 
by the mean-field theory, which is definitely a better approximation at scales of the order 
of the final gap. The one-loop flow without order parameter feedback fails even in pure 
mean-field models at the scale where symmetry breaking sets in. Hence we have to stop 
the one-loop flow at a scale safely above the scale of symmetry breaking. On the other 
hand one would like to capture the one-loop fluctuations not treated in mean-field theory 
down to the lowest possible scales. Unfortunately there is no unique choice of an "optimal" 
Amf. In the following we will stop the one- loop flow at a scale Amf which is twice as big 
as the largest (as a function of k) gap value in the symmetry-broken state obtained from 
the mean-field theory for the states below Amf- The relation Amf = corresponds to 
the straight line in Fig. 8. The estimate for is obtained from the intersection of that 
line with the largest Ak, that is, the ones with k on patch 1. In the present case this 
estimate is rather robust with respect to slight shifts of Amf- 

B. d-wave superconductivity 

It is well known that the one-loop RG flow generates an attractive interaction in 
the d-wave Cooper channel.— 'i^'^^'i^ Although antiferromagnetic fluctuations contribute 
signiflcantly to this attraction, the Cooper instability ultimately competes with antifer- 
romagnetism. In Fig. 9 we show results for the superconducting gap function Ak as a 
function of the cutoff Amf in the RG+MF scheme. Different curves correspond to different 
patches, as in Fig. 8. The parameters have been chosen such that the antiferromagnetic 
interaction is too weak to drive antiferromagnetic order. Within plain mean-fleld the- 
ory no superconductivity is obtained since the bare (repulsive) Hubbard interaction is 
repulsive in the s-wave Cooper channel and vanishes for any other symmetry. Hence Ak 
vanishes for Amf at the band edge and increases monotonically upon lowering Amf- The 
pronounced momentum dependence of Ak reflects the d^'^-y'^ symmetry. As in the case 
of antiferromagnetism, Ak diverges if the one-loop flow is continued below scales of the 
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FIG. 9: (Color online) Superconducting g function of the cutoff Amf for the Hubbard 

model in a regime where superconductivity is the only instability. Parameters: U = 2.5, t' = 
-0.2, n = -0.9265. 

order of the gap amplitude. A reasonable choice of Amf is obtained from the condition 
Amf = 2Ak (straight line in Fig. 9) applied to the largest Ak (patch 1). As long as one 
does not approach the divergence, the result for Ak does not depend too much on the 
precise choice of Amf- 

C. Antiferromagnetism and superconductivity 

We now analyze the competition between antiferromagnetism and superconductivity 
within the full RG+MF theory, where both types of order and also the possibility of their 
coexistence are allowed. We focus on the case t' < 0, which is realized in particular in the 
cuprate superconductors. We restrict ourselves to densities at and below half-filling. 

In Fig. 10 we show results from the RG-I-MF calculation for the amplitudes of the 
d-wave superconducting gap Ak and the antiferromagnetic gap Ak as a function of the 
chemical potential /i, and in Fig. 11 as a function of the electron density. The latter is 
computed simultaneously with the order parameters and differs from the bare density 
corresponding to fi. The interaction U = 2.5 is strong enough to stabilize an antifer- 
romagnetic insulator at half-filling, in spite of the magnetic frustration induced by the 
next-to-nearest neighbor hopping t' = —0.15. The system is fully gapped at half-filling 
and the superconducting order parameter is strictly zero there. With decreasing filling 
the antiferromagnetic gap decreases monotonically, while Ak remains numerically zero. 
For electron densities below one, holes appear first in pockets around (7r/2,7r/2), which 
define a surface of low-energy excitations of the non-half filled system. In principle, this 
residual Fermi surface is always unstable against superconductivity, due to the attrac- 
tive interaction in the d-wave Cooper channel. However, very close to half-filling, where 
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FIG. 10: (Color online) Amplitudes of the superconducting gap and the antiferromagnetic 
gap as a function of fx for U = 2.5 and t' = —0.15. The dotted lines are the separatrices 
between different Fermi surface topologies as specified in Fig. 6. 
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FIG. 11: (Color online) Amplitudes of and as a function of density for U = 2.5 and 
t' = —0.15. Filled colored symbols represent the results obtained for the combined theory with 
two order parameters, while in the results represented by open symbols only one order parameter, 
either antiferromagnetic or superconducting, was allowed in the mean-field calculation. The two 
slightly different results for the antiferromagnetic gap at half-filling correspond to two different 
choices of fi within the gap; this /^-dependence is an artefact of the approximations. 



the pockets are small, the superconducting gap is tiny, since the d-wave attraction is 
very small near the Brillouin zone diagonal. When the hole pockets are large enough, 
roughly for n < 0.92, a sizeable superconducting gap develops, coexisting with a finite 
antiferromagnetic order parameter. When fi approaches /IvH; the antiferromagnetic gap 
decreases rapidly. This leads also to a very fast decrease of the electron density (see 
Fig. 11) as a function of decreasing /i. The antiferromagnetic gap vanishes for /x < /x^h, 
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while Ak remains finite and decreases monotonically with decreasing filhng. Note that 
the antiferromagnetic transition is not generically situated as close to /ivH as in Fig. 10. 

The numerical analysis is rather involved in the coexistence region and we have not 
clarified the nature of the magnetic transition (first or steep second order). The presence 
of the superconducting gap makes the magnetic transition smoother, but whether this 
effect is big enough to yield a continuous transition is not clear from the numerical data. 
In any case, fiuctuations which have been neglected in the mean-field approximation are 
expected to play an important role in the transition region and may even affect the order 
of the transition. 

The respective results for each order parameter when the other one is set to zero are 
also shown in Fig. 11. When is set to zero, the superconducting gap Ak persists 
even at half filling. When A^ is set to zero, the antiferromagnetic order parameter is 
enhanced in the coexistence region. In both finite value for one order parameter 

leads to a suppression of the other. While these results have been obtained for a weak 
on-site repulsion, they are in line with the behavior at stronger coupling as obtained 
by the variational Monte Carlo technique,— by cluster and cellular dynamical mean- 
field theory,^-'*^- and by variational cluster approximations^"^ for the two-dimensional 
Hubbard model. The regime of substantial coexistence of antiferromagnetic order and 
superconductivity obtained in those calculations is typically larger than in our results. 
However, in the cluster calculations the amount of coexistence seems to depend sensitively 
on the cluster size.— In a recent cellular dynamical mean-field calculation with a 2x2 
cluster it was found that coexistence is not possible if U exceeds the band width.— 

In Fig. 12 we show the angular dependence of Ak and for three selected densities, 
along with the corresponding effective Fermi surfaces. The angle is defined with respect 
to the kx axis. For n = 0.99 we obtain a sizable antiferromagnetic gap which has s- 
wave symmetry and is slightly anisotropic. The superconducting gap is practically zero. 
The hole pocket enclosed by the effective Fermi surface is rather small and does not 
support a sizable superconducting gap. For n = 0.906 we observe a coexistence of both 
order parameters. While Ak is reduced in size compared to the case n = 0.99, its shape 
remains essentially the same. The hole pocket extends further away from the Brillouin 
zone diagonal and allows for a substantial superconducting gap with d-wave symmetry. 
In this situation the low-energy excitations are gapped due to Ak near the points (vr, 0) 
and (0,7r), and due to Ak along the hole pocket, with nodes along the Brillouin zone 
diagonal. For the even smaller density n = 0.862 the antiferromagnetic order parameter 
vanishes and only a d-wave superconducting gap remains, which extends over the whole 
Fermi surface except at the nodal points. Note that the momentum dependence of the 
gap function has peaks near the van Hove points and is relatively flat near the Brillouin 
zone diagonal, implying that terms beyond the lowest d-wave harmonic COS kx cos Jvy 
contribute significantly to Ak- 
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k 

X 



FIG. 12: (Color online) Top: Angular dependence of the d-wave superconducting and antifer- 
romagnetic order parameters at three different values of the chemical potential. Bottom: Fermi 
surfaces in the magnetic Brillouin zone (solid blue lines); in the antiferromagnetic state close 
to half-filling the Fermi surface forms a hole pocket around (7r/2,7r/2). The corresponding bare 
Fermi surfaces, backfolded with respect to Umklapp surface, are shown as broken lines. The 
straight lines indicate the patching scheme. The parameters U = 2.5 and t' = —0.15 are the 
same as in Figs. 10 and 11. 

In Figs. 10-12 the coupling strength U has been chosen such that three different states, 
that is, antiferromagnetic insulator, doped antiferromagnet coexisting with superconduc- 
tivity, and a pure d-wave superconductor could be obtained for different choices of fi. By 
contrast, for a sufficiently small U (at finite t') the system is a d-wave superconductor for 
any /i. On the other hand, for large U the system is an antiferromagnet near half-filling 
and switches to a purely superconducting state if the chemical potential descends below 
a certain critical value. The superconducting gap in the doped antiferromagnet remains 
tiny for large U. A global view of the U and /i-dependences is presented in Fig. 13, where 
we show the ground state phase diagram of the two-dimensional Hubbard model in the 
fi-U plane for a fixed t' = —0.1. Note that the parameter region supporting substantial 
coexistence of antiferromagnetism and superconductivity is rather narrow. For U = 1.75 
and U = 2, antiferromagnetic order is stabilized only below half-filling. 



23 



2.5 r 



U 



1.5 



1^ 



o 


sc 


o 


AF 




coex 



-0.4 



-0.3 



-0.2 



FIG. 13: (Color online) Points with antiferromagnetism, superconductivity and coexistence of 
both in the fj,-U plane for a fixed t' = —0.1; "coexistence" means that the amplitude of Ak/A^ 
is at least 10^'^. 



V. CONCLUSION 



We have analyzed the competition between antiferromagnetism and superconductiv- 
ity in the two-dimensional Hubbard model by combining a functional RG flow with a 
mean-field theory for symmetry-breaking. Effective interactions were computed by in- 
tegrating out states with an energy above a scale Amf via a one-loop RG flow, which 
captures in particular the generation of an attractive interaction in the d-wave Cooper 
channel from fluctuations in the particle-hole channel.— >^>^ These effective interactions 
were then used as an input for a mean-field treatment of the remaining states, below the 
scale Amf, with commensurate antiferromagnetism, singlet superconductivity, and triplet 
TT-pairing as the possible order parameters. Triplet vr-pairing appears generically when 
antiferromagnetism and singlet superconductivity coexist.^'* Our theory extends previous 
mean-field treatments of antiferromagnetism and d-wave superconductivity, where the 
effective interactions were specified by a relatively simple ansatz with a few input pa- 
rameters, or they were identified with bare interactions of a microscopic model.— "^"^"^ 
In the present work the size and shape (momentum dependence) of the effective interac- 
tions were actually computed, within an approximation that captures fluctuations and is 
controlled at weak coupling. 

It turned out that the feedback of vr-pairing on the other order parameters is negligble 
for the two-dimensional Hubbard model. The key players, antiferromagnetism and d- 
wave superconductivity, strongly compete. For a sufficiently large U (depending on the 
size of t') an antiferromagnetic insulator is stabilized at half-filling, as expected. Doping 
the antiferromagnet with holes (for t' < 0) leads to an effective Fermi surface with hole 
pockets around [it /2,tt/2). In principle, these give rise to a superconducting instability. 
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but the corresponding gap is usually tiny. At larger doping the antiferromagnetism breaks 
down and d-wave superconductivity prevails. There is a small range of densities where 
both orders can coexist with a sizable order parameter for each. 

Coexistence of antiferromagnetism and superconductivity in the two-dimensional 
Hubbard model has also been obtained by the variational Monte Carlo technique,— 
by cluster and cellular dynamical mean-field theory,—*^ and by variational cluster 
approximations.—"^ Coexistence has also been observed in some cuprate high-temperature 
superconductors in the underdoped regime.— 

We have restricted our analysis of symmetry breaking to superconductivity and com- 
mensurate antiferromagnetism, while in principle also other instabilities may arise. The 
most serious candidates (and difficult to treat) are probably the various possibilities of 
incommensurate magnetic order in the low doping regime.— For the two-dimensional 
t-J model with a finite t' it was shown that for very low doping a spiral phase is stabilized 
and that superconductivity develops in the spiral state.'^ Another possibility is the for- 
mation of a d-density wave, that is, a charge density wave with a wave vector near (vr, tt) 
and a form factor with d-wave symmetry.— This ordered phase is known to occur in a 
suitable large extension of the t-J model, where it competes and partially coexists with 
d-wave superconductivity.— Renormalization group calculations for the two-dimensional 
Hubbard model yield enhanced d-density wave correlations, but the antiferromagnetic 
and superconducting instabilities are stronger.§S 

A more complete analysis of the two-dimensional Hubbard model by renormalization 
group methods, including also order parameter fiuctuations at low energy scales, can be 
expected to yield further important clues for a better understanding of the interplay of 
magnetism and superconductivity. 



We thank R. Gersch, C. Honerkamp, A. Katanin, M. Salmhofer, P. Strack, O. Sushkov, 
and R. Zeyher for valuable discussions. 



Let 0k be any of the expectation values (riko-), (^ko-)' (Pk), (Pk), or their complex 
conjugates. Applying the necessary condition for a minimum, dfl/d(j)]i = 0, to f2 in Eq. 
(EOD yields 
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APPENDIX A: DERIVATION OF GAP EQUATIONS 
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The derivatives of the energy eigenvalues E^, can be expressed in terms of minors 
det(M^;-'-' ) of M^, and matrix elements (A^k')ij' of -^k' as 

50. -E.det(M,y)f^ 50. ( ^^H^k^ )■ (A2) 

This follows from acting with d/ 90. on the equation for the eigenvalues det (^Aik'—E^, l) = 
and the rule for derivatives of determinants 

^det(M) = J2 (-1)^"^' det(M^-^') . (A3) 



The gap equations are now obtained by taking the specific derivatives for each case. 
Terms differing only by a momentum shift Q can be summed up by using the identities 
E^^Q = and det(M^:j?'Q = det(Mk'-''+^'^''+^) , with j + 2 and f + 2 modulo 4, which 
follow from the block structure ( |2711 of A^k- 

To derive Eq. (136|) we start from Kramer's rule for matrix inversion in terms of 
minors, det(M) (M^^)ii' = (— l)'~''det(M'''). Summing the diagonal elements yields 

det(M'') = det(M) tr(M^^). Expressing the trace and the determinant by the eigen- 
values rua of M yields E;det(M") = J^aYla'T^a^a' ■ If one of the eigenvalues, say mx, 
vanishes, only the term a = X contributes, yielding E;det(M") = Yly^/^x^y- Applyi^g 
this identity to the matrix one obtains directly Eq. fl36|) . 



APPENDIX B: ATTRACTIVE HUBBARD MODEL 



In the attractive Hubbard model away from half-filling only the Cooper channel leads 
to an instability, namely s-wave superconductivity.™ All other channels yield just finite 
renormalizations. In Fig. 14 we show the dependence of the superconducting gap on Amf 
as obtained from the RG+MF theory for the attractive Hubbard model with U = —1.5 
and t' = —0.1 at quarter- filling. As in the case of antiferromagnetism in the repulsive 
model, the gap size is reduced by fluctuations compared to the mean-field result. For A 
between 4 — f» 6 and A + n ^ 2 only states above the Fermi level are integrated out, while 
states below the Fermi level are captured only for smaller A. This is the reason for the 
kink in at Amf = 4 + /i. The gap diverges at a critical scale A^ due to the breakdown 
of the one-loop approximation discussed already in Sec. IV A. For Ac <^ Amf 1 the 
dependence of Ak on Amf is relatively weak. This is the regime where the cutoff is still 
high enough so that the one-loop approximation is still accurate, but on the other hand 
already sufficiently low so that the flow of the effective interaction Vkk' is dominated by the 
particle-particle channel, such that interactions beyond the reduced (mean-field) model 
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FIG. 14: (Color online) Superconducting gap Ak as a function of the cutoff Amf for the at- 
tractive Hubbard model {U = —1.5, t' = —0.1) at quarter-filling. Due to the weak momentum 
dependence of curves corresponding to different patches lie almost on top of each other. The 
inset shows A^ for Amf ^ 1 with a higher resolution. 

are irrelevant. 
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